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ABSTRACT 

High energy observations of extreme BL Lac objects, such as 1ES 0229+200 or 1ES 0347- 
121, recently focused interest both for blazar and jet physics and for the implication on the 
extragalactic background light and intergalactic magnetic field estimate. However, the number 
of these extreme highly peaked BL Lac objects (EHBL) is still rather small. Aiming at increase 
their number, we selected a group of EHBL candidates starting from the BL Lac sample of 
iPlotkin et al.l(l201 ll) . considering those undetected (or only barely detected) by the Large Area 
Telescope onboard Fermi and characterized by a high X-ray vs. radio flux ratio. We assem¬ 
bled the multi-wavelength spectral energy distribution of the resulting 9 sources, profiting of 
publicly available archival observations performed by Swift , Galex, and Fermi satellites, con¬ 
firming their nature. Through a simple one-zone synchrotron self-Compton model we estimate 
the expected VHE flux, finding that in the majority of cases it is within the reach of present 
generation of Cherenkov arrays or of the forthcoming Cherenkov Telescope Array (CTA). 

Key words: galaxies: active - galaxies: jets - radiation mechanisms: non-thermal - gamma- 
rays: galaxies. 


1 INTRODUCTION 

Intense emission of 7 rays is a distinctive feature of blazars, ac¬ 
tive galactic nuclei (AGN) dominated by the boosted non-thermal 
continuum from a relativistic jet pointed toward the observer. The 
2LAC catalogue (Ackermann etal, 201 lj), listing the AGN detected 
with high significance by the Large Area Telescope onboard Fermi 
during its first two years of operations, contains 886 sources, of 
which 862 are blazars (395 sources classified as BL Lacs, 310 Flat 
Spectrum Radio Quasars (FSRQ), and 157 sources of “unknown 
type”). Blazars also dominate the extragalactic sky at very high 
energy (VHE, E > 100 GeV) with BL Lac being the dominant 
population (58 over a total of 67 extragalactic sources discovered 
until November 2014 according to TeVCaQ). Within BL Lacs, the 
large majority (51 out of 58) of the VHE emitters belongs to the 
High-peaked BL Lac (HBL) subclass. 

The spectral energy distribution (SED) of blazars displays 
two characteristics broad humps, who se peak frequency appears 
to anticorrelate with the emitted power dFossati et al Jl 1998L but see 
iGiommi et al.ll2005h . While the low energy (from radio up to op¬ 
tical) emission is clearly associated to synchrotron radiation, the 
nature of the mechanisms responsible for the high-energy contin¬ 
uum is still debated. The majority of the studies adopts the so called 
leptonic scenario, in which the high energy radiation is interpreted 
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as the product of the inverse Compton (IC) scattering of the rel¬ 
ativistic electrons and soft photons (either produced internally, i.e. 
the synchrotron photon themselves, or externally to the emitting re¬ 
gion). Hadronic models dBottcher et alJ[2ol3l) instead assume that 
the 7-ray emission is either the byproduct of reactions initiat ed by 
ultrarelativistic hadrons and m ediated by mesons (e.g. lMticke et al.1 
l2003l : lAtovan&Den]igj 2d03l) or direct synchrotron emission from 
protons (e.g. Aharonianll200d) . 


For the great majority of known BL Lac objects the high- 
energy component peaks in the 1-100 GeV band. However, 
there is a small group of sources detected at TeV energies by 
Cherenkov arrays for which the maximum is located above 1 
TeV (e.g. iTavecchio et al.l 1201 11) . The extreme hardness of the 
spectrum makes these sources very faint and thus often unde¬ 
tected in the Large Area Telescope (LAT) band. These features 
can be exploited to effectively constrain the intergalactic magnetic 


field (IGMF, e.g. Neronov & Vovk 2010; Tavecchio et al 

2010 

20l]l; Dolag et al. 2009; DermeretalJ 201 l|T Vovk et al. 

2012 

Oikonomou et al. 2014|). Other properties shared by these sources 


are the extremely large ratio between the X-ray and the radio flux 
and the hardness of the X-ray continuum (F ~ 2), locating the syn¬ 
chrotron peak in the medium-hard X-ray band. These characteris- 
tics lead to collect them under the term “extreme” HBL (EHBL, 
ICostamante et al.ll200ll) . In the IR-optical regime the emission is 
dominated by the host galaxy and the non-thermal jet continuum 
starts to be important only in the UV band. The modeling of their 
SED within the Synchrotron Self-Compton (SSC) framework re- 
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Figure 1. Spectral en ergy distribution of 1ES 0229+200 (left) and 1ES 0347-121 (right), two of the most representative EHBL det ected at TeV energ ies (see 
iTavecchio et all20 11] for references). Blue symbols show the TeV spectrum corre cted for the absorp tion by EBL using the model of lDormnguez et a l] J2011|). 
The black points for 1ES 0229+200 report the Fermi/LAT spectrum obtained bv lYovk et alj feOlit) , while those for 1ES 0347-121 come from lTanaka et alj 

iioll . 


veals rather unusually l ow magnetic fields (B < 0 .01 G) and 
large electron energies ( ITavecchio et al.l 1201 Ot 1201 lh . The emit¬ 
ted hard TeV spectrum (once corrected for interaction with the 
extragalactic background light, EBL) can be reproduced assum¬ 
ing that the electron energy distribution is truncated below a mini¬ 
mum Lorentz factor around 7 m j n ~ 1 0 5 dKatarzvnski et al .1120051 : 
Tavecchio et ai1l2009l . 1201 it iKaufmann et al.ll201lh . This interpre- 
tation also accounts for the peculiar UV/X-ray spectrum. Alter¬ 
natively, the SEP could be the result of a Max wellian electron 
distri b ution iLefa et al.l20 lH), inte rnal absorption dAharonian et al.1 
120081 : IZacharopoulou et afll BOl lh or inverse Compton scattering 
between electrons and photons of the cos mic microwave back¬ 
ground in the large-scale (~kpc) jet dBottcher et al.ll2008l) . Alter¬ 
natively lepto-hadronic models can be also invoked, either through 
proton-synchrotron emission, or through secondary cascades inside 
the emi ssion region initiated by ultra-relativistic hadrons (see e.g. 
ICerruti et alj2015l ). A last suggestive possibility is that high-energy 
photons are produced by ultra-high energy protons along the line of 
sight injected by the blazars into the intergalactic space 1 Essev et al.l 

1201 ll : iMurase et al.ll2012al : Izheng & Kandl20l3klTavecchioll2014 ). 

From the brief description given above is clear that EHBL are 
rather interesting objects, both for the study of jets phenomenol¬ 
ogy or even Ultra High Energy Cosmic Ray (UHECR) astrophysics 
and for the use of probes of EBL and IGMF. However, their use is 
somewhat hampered by the small number of EHBL detected at TeV 
energies. This is the main fact driving the present work, in which 
we intend to characterize a group of EHBL detectable at TeV en¬ 
ergies by present instruments or by the upcoming Cherenkov Tele¬ 
scope Array (CTA. lAcharva et aHl2013l) . Sp arse groups of EHBL 
have been identified in previous work (e.g. iGiommi et alj 1200a : 
iNieppola et al.1 12OO6I) . Here we intend to follow a focused and 


well defined selection procedure, based on the compilation of the 
SDSS/FIRST BL Lac of [piotkin et al.l j201 ll) . and the criterion of 
an extreme radio-to-X-ray flux ratio (§2). We show their SED using 
recent Swift data (§3). We also profit from archival Galaxy Evolu¬ 
tion Explorer (GALEX) satellite data and we analyze the whole 
FermilLAI photon archive in order to constrain their behavior in 
gamma rays by computing at least upper limits if no detection is 
possible, which comes out to be the most common result, in agree¬ 
ment with our expectations. Based on a simple synchrotron self- 
Compton homogeneous model we reproduce the resulting SED 
(§4), giving also an estimate of their flux and detectability for Imag¬ 
ing Atmospheric Cherenkov Telescopes (IACT) of the present and 
forthcoming generation. 

Throughout the paper, we assume the following cosmological 
parameters: Ho = 70 km s _1 Mpc~\ Bm = 0.3, Ha = 0.7. We 
use the notation Q = Qx 10 A in cgs units. 


2 SELECTION OF TEV CANDIDATES EHBL 

The SED of two representative EHBL detected in the TeV band, 
1ES 0229+200 and 1ES 0347-121, are reported in Fig. [fl To se¬ 
lect EHBL candidates we are guided by two evident peculiarities 
of these SED, namely the large X-ray/radio flux ratio and a non- 
thermal optical continuum lower than the thermal contribution from 
the galaxy. These two features suggest to select EHBL among the 


2 Note that due to a bug in th e calcu lation of the effective area, the X- 
ray spectrum reported in Tavecch io et_al.j 120091 ) was too high. In Fig.[7]the 
correct spectrum is shown. 
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BL Lacs with high X-ray/radio flux ratio whose optical spectrum is 
dominated by the galaxy. 

We emphasiz e that this c riterion should not be confused with 
the one adopted in lCostamante & Ghisellinil ( i2002l) who aimed to 
maximize the TeV flux (with a selection based on the evidence 
of large X-ray and radio fluxes), while here we tend, in a wide 
sense, to minimize the TeV spectral index. The criterion is rather 
model independent, as it arises by similarity with the archetypal 
EHBL 1ES 0229+200 and 1ES 0347-121 mentioned above. In the 
scope of a pure one-zone SSC model (e.g.. lTavecchio et~al]| 19981) . 
this can be interpreted as the outcome of an electron distribution 
characterized by high minimum ( 7 min) and break ( 7 b) Lorentz fac¬ 
tors; this translates into large values of the peak frequencies for 
the synchrotron (v s ) and Inverse Compton (vie) components of 
the radiation spectrum, and dim ( due to scarc ity of seed photons) 
but hard TeV spectra. Icostamante & Ghisellinil I 2 OO 2 I) instead re¬ 
quested both high energy electrons and seed photons to be abun¬ 
dant in the emission region, in order to maximize the bolometric 
VHE luminosity. This interpretation suggests that the regions of the 
Fx /F r parameter space plotted in Fig.[2]selected by the two crite¬ 
ria should actually host BL Lac populations of different flavour, 
with the TeV-brightest clustered towards the upper-right corner 
while the Te V-hardest clu stered at the upper-left corner. While the 
ICostamante & Ghisellinil i2002h selection has been robustly proved 
by the detection at TeV energies of most of their candidate sources, 
the one we adopt here will need extensive observation of the se¬ 
lected sources in order to be proven or rejected. 

Guided by this criterion, we started our selection from the list 
of 71 BL Lac presented bv lPlotkin et alj <1201 lh resulting from the 
correlation of SDSS and FIRST surveys and by optical spectrum 
dominated by the host galaxy emission; to enforce this, we limit 
the study within z = 0.4. To apply our criteria we further selec t 
those with measured X-ray flux (by ROSAT), from IPlotkin et alj 
d2010h . obtaining a total of 50 BL Lacs. Note that the cut in redshift, 
originally dictated by the requirement of the galaxy dominance in 
the optical band, also ensures that the relatively small suppression 
of the EBL for 7 rays up to a few TeV does not prevent detection 
by current Cherenkov arrays. 

The X-ray flux in the 0.1-2.4 keV band and the radio flux 
at 1.4 GHz (expressed as Fr = u¥ v for homogeneity with the 
integrated X-ray flux), of the sources resulting from the selec¬ 
tion have been taken from iPlotkin et alj d2010h and are reported 
in Fig|2] There, the blue points show the sources already detected 
by LAT (2FGL catalogue). For comparison, the magenta points re¬ 
port the fluxes of 1ES 0229+200 and 1ES 0347-121 (although not 
belonging to our selected sources). The oblique dashed line corre¬ 
sponds to a fixed ratio Fx/Fr = 10 4 . This particular, rounded 
but somewhat arbitrary value is suggested by the position of the 
two known EHBL in the plane. The sources above the line are ex¬ 
pected to be good EHBL candidates, without any claim for com¬ 
pleteness, as we are not aiming at it at this level. For instance, 
just outside of the region and next to 1ES 0347-121 lies RBS 
1049, which would easily fit in the same scheme. In this region 
there are two sources already detected by LAT, RBS 0723 and 
1ES 1426+428. The latter is already a well established TeV source 
dHoran et alj2Qo3) . Although its SED is coherent with our picture, 
it can be considered a transition object towards the “bright X-ray”- 
“bright radi o” c orner of the plot, populated by the bright blazars 
of lCostamante & Ghisellinil d2002l) therefore will be omitted in the 
following discussion. It’s worth noticing that 1ES1426+428 is not 
a standalone case, as another hard s ource not belonging to our sam¬ 
ple but to the one of lGiommi et alj d2Q05l ) that is built along a sim¬ 


ilar line has been discovered atWHE by HESS (SHBL J001355.9- 
185406 at 2 = 0.095. IH.E.S.S. Collaboration et alJlioIj) . More 
recently RBS 0723 has been discovered at VHE too, by MAGIC 
dMirzovan 2014 ), Other extr eme BL Lacs such as those contem¬ 
plated in Cerruti et al.1 d2015l) also show X-ray to radio flux ratios 
high enough to fit in the scheme. 

In Tab. [I] we report the names of the selected EHBL together 
with their coordinates and redshifts. Note that, as an indirect result 
of the selection, requiring a bright X-ray emis sion, most of them 
belong to the Rosat Bright Survey (RBS. fvbges et alJI 19991) . 


3 MULTI-WAVELENGTH DATA 

For the sources fulfilling our criterion on the Fx /Fr ratio we built 
multiwavelength SED exploiting publicly available data. We used 
archival optical-UV data from Swift !UVOT and soft X-ray data 
from Swift/X RT for the vast majority of the sources, while for RBS 
0921 and RBS 1176 we requested dedicated Swift observations. We 
also profited from the GALEX archivfl Other data have been also 
added using the ASI Science Data Center (ASDC) archivtQ- All 
the considered sources have been observed at least one time with 
Swift, ensuring a good description of the crucial UV-X-ray band. 
No particular attempt was made to seek for strictly simultaneous 
observations, that would need a dedicated observational campaign, 
being these sources generally reputed of secondary importance and 
therefore left aside from any monitoring campaign of some time 
density. 

We also analyzed FermitLKF data. Almost all the sources are 
undetected and we could only calculate upper limits. 

In the following we describe the analysis performed on the 
Swift/X RT Swift !UVOT and Fermi! LAT data. 


3.1 Swift/XRT data 

Swift/X RT ^Burrows et al.l200l) observations were available for all 
the sources in our subsample, except for RBS 0923 and RBS 1176, 
which were targets of dedicated observations: RBS 0923 was ob¬ 
served in July 2012, while RBS 1176 was observed in November 
2011, then in May and June 2012 . For most of the other sources 
only one observation was available, except for RBS 1029 which 
was observed in November and December 2007, 1ES0927+500, 
observed in September 2010 and March 2011, and RBS 1510 in 
June, September, October 2011 and in January 2012. 

Swift/XRT data were analyzed by using the HEASoft v. 6.13 
software package with the CALDB updated on 21 January 2013 
and processed with xrtpipeline v. 0.12.6 with standard 
parameters. Spectra have been grouped to have at least 20 counts 
per bin, in order to use the \ 2 test and analyzed with xspec v. 12 . 8.0 
in the 0.3-10 keV energy band. 

For all sources (with perhaps the exception of RBS 1176, see 
below) an absorbed power law model provides a good description 
of the spectrum. In the majorit y of ca ses the ab sorp tion column can 
be fixed to the Galactic value jKalberla et al.ll2005h . 

In Tab. [2] we report the results of the fitting procedure. In the 
two cases of 1ES 0927+500 and RBS 1029 we separately consider 
two spectra. The best fit parameters are perfectly consistent. All 


3 http://galex.stsci.edu/galexview/ 

4 http://tools.asdc.asi.it/ 



















































4 Bonnoli et al. 


i 

00 

CV2 

I 

£ 

o 

GO 

<D 


> 

0 

44 

cvi 


d 


X 

Pxh 

GO 

o 

hJ 


-10 


-11 


-12 


-13 


- 

1 ' ' ' ' 1 

1 1 1 1 1 1 1 

1 1 1 1 1 1 ! 

1ES 1426+429.'' 

• 


- 


1ES 0927+500 


“ 



O 0723 


“ 


• 


“ 



1ES 0229+200 




0921 

- 1510 




O 


• 




1029 •••' 

o 

c 


- 

15 <g o 0923 O 

• oo° 

o 

- 


1176 

O 



- 


• • o 


- 



1ES 0347-121 



_BZB 

J0022 + 0006 

c 9 

° o 


O 


o o 




m oc 
o o 
o 

o 

o 

c 

o 

• 

o o 



O 

o 


o 

o 

- 

- 


o 

o 

- 

- 



• 

- 

- 

1 ■ ■ ■ ■ 1 

■ ■ ■ ■ 1 ■ ■ 

■ ■ 1 ■ ■ ■ ■ 



-16 - 15.5 -15 - 14.5 -14 

Log u r F r t er S cm ' 2 s_1 ] 


Figure 2. Fx vs. Fr diagnostic plot for the sample of BL Lacs described in the text. Red open circles represent the full sample; the blue filled circles are the 
LAT detected ones. The 12 sources above the black dotted diagonal line have a high ( Fx / Fr > 1 0 4 ) ratio of X-ray vs. radi o flux. T he magen ta filled circles 
are the two archetypal TeV detected, but GeV faint extreme BL Lacs, 1ES 0229+200 tVovk et all20121) and 1ES 0347-121 ^Tanaka et alJl2014h . 


Source Name 

R.A.(J2000) 

<5(J2000) 

1 

b 

Redshift 

A B 

BZB J0022+0006 

5.5040 

0.1161 

107.18 

-61.85 

0.306 

0.108 

RBS 0723 

131.8039 

11.5640 

215.46 

30.89 

0.198 

0.093 

1ES 0927+500 

142.6566 

49.8404 

168.19 

45.71 

0.187 

0.073 

RBS 0921 

164.0275 

2.8704 

249.28 

53.28 

0.236 

0.178 

RBS 0923 

164.3462 

23.0552 

215.96 

63.91 

0.378 

0.088 

RBS 1029 

176.3963 

-3.6671 

273.11 

55.34 

0.168 

0.130 

RBS 1176 

193.2540 

38.4405 

121.36 

78.68 

0.371 

0.083 

RBS 1510 

233.2969 

18.9081 

29.21 

52.05 

0.307 

0.210 

RBS 1555 

241.3293 

54.3500 

84.35 

45.60 

0.212 

0.041 


Table 1. List of the extreme blazar candidates selected from the sample of fpiotkin et al] | 2011|) . For e ach source the equatorial (J2000) and galactic coordinates 
are reported (in degrees), the redshift and the Ag extinction coefficient from fschlegel et al ]| 1998h . As a consequence of the selection within SDSS, all the 
sources lye at high galactic latitudes. 


the other multiple datasets could be merged as no hint of variability 
was found. 

More complex is the case of RBS 1176. A fit of the summed 
spectrum with a power law model returns an absorption column 
Nh largely in excess to the Galactic value (Table 2). Alternatively, 


a good fit ( \ 2 /d.o.f. = 12.46/13) can be obtained assuming a 
broken power law model and the Galactic value of Nh, with slopes 
IT = 1.1 ± 0.3, T 2 = 2.3 ± 0.3 and break energy Eb = 1.3 ± 0.3 
keV. The spectra of the single pointings are much less constrain¬ 
ing. A single power law model (with photon index T ~ 1.7) with 
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Source Name 

Obs. ID 

Exp. Time 
[s] 

r 

N h 

[10 20 cm -2 ] 

X 2 /d.o.f. 

Fo.3-10 keV 
[10 —12 erg cm -2 s—1] 

BZB J0022+0006 

38113001 

4700 

2.40 ±0.25 

2.76 (Gal.) 

7.19/8 

1.7 ±0.3 

RBS 0723 

37396001 

2000 

1.78 ±0.12 

3.17 (Gal.) 

20.8/18 

2.6 ±0.3 

1ES 0927+500 1 

39154001 

2870 

2.0 ±0.1 

1.38 (Gal.) 

40.13/37 

14 ±1 

2+3 

39154002-3 

2260 

2.2 ±0.1 

1.38 (Gal.) 

29.4/33 

16 ±1 

RBS 0921 

37547001 

4700 

1.89 ±0.08 

3.82 (Gal.) 

44.8/46 

8.4 ±0.4 

RBS 0923 

48001001-4 

6400 

2.2 ±0.1 

1.12 (Gal.) 

25.8/26 

3.6 ±0.3 

RBS 1029 1 

36813001 

2900 

2.2 ±0.2 

2.22 (Gal.) 

8.1/6 

2.7 ±0.3 

2 

36813002 

2900 

2.0 ±0.1 

2.22 (Gal.) 

9.35/10 

3.2 ±0.2 

RBS 1176 

48000001-3 

5250 

2.3 ±0.3 

13.3 ±6.5 1 

14.1/13 

2.7 ±0.3 

RBS 1510 

91101001-5 

3700 

2.2 ±0.1 

3.83 (Gal.) 

35.0/33 

8.3 ±0.6 

RBS 1555 

38303001 

7100 

2.0 ±0.1 

0.886 (Gal.) 

23.4/29 

1.7 ±0.1 


Table 2. Results of the Swift/XKT data analysis for all the archival observations of our sample available until the end of August 2012. The table reports: source 
name, exposure time (expressed in s), best fit value of spectral index of power-law model, column density (in units of 10 20 cm -2 ), y 2 and degrees of freedom 
of the best fit, flux in the 0.3-10 keV band (in units of 10 -12 erg cm -2 s -1 ). L For this source N Hga i = 1.72 X 10 20 cm -2 . 


Galactic absorption is barely compatible with the data but a curved 
spectrum is clearly suggested by the shape of the residuals. 

3.2 Swift/WOT data 

Swift /UVOT l lRoming et al.ll2005l ) is a 30 cm diffraction-limited 
optical-UV telescope, equipped with six different filters, sensitive 
in the 1700-6500 A wavelength range, in a 17’ x 17’ FoV. We 
retrieved from the HEAS ARC database the UVOT images in which 
our target sources were observed. The maximum angular distance 
from the optical axis does not exceed 4’ for any source. 

For all the sources and the available different observations the 
analysis was performed with the uvotimsum and uvotsource 
tasks with a source region of 5". while the background was ex¬ 
tracted from a source-free circular region with radius equal to 50”. 
The extracted magnitudes were corrected for Galactic extinction 
using the values of Schleg eTet al.l i ll9981) (repo rted in the last col¬ 
umn of Tab. |T] and applying the formulae by IPeil 1 1 9921) for the 
UV filte rs, and eventually were converted into fluxes following 
IPoole et alj fcOOSh . 

Tab. [3] reports the observed Vega magnitudes in the 
Swift/UVOT v, b, u, ml, m2, and w2 filters, together with statis¬ 
tical uncertainties. Systematic uncertainties are never greater than 
0.03 mag and therefore dominated by statistical ones in the vast 
majority of cases. 

3.3 Fermil LAT data 

Publicly available FemiifLAT data were retrieved from the Fermi 
Science Support Center (FSSC) and analyzed by means of the LAT 
Science Tools v. 9.27.1, together with the Instrument Response 
Function (IRF) Pass 7 and the corresponding isotropic and Galac¬ 
tic diffuse background models. Source (class 2) photons in the 0.1- 
100 GeV energy range, collected until 21 August 2012 and com¬ 
ing from direction within 10 degrees from the nominal position 
of the source were selected and filtered through standard FSSC 
quality cuts. Standard analysis steps were then performed, even¬ 
tually adopting the test statistic from IMattox et al.l i 1996b to assess 
the significance of excess signal in correspondance with our tar¬ 
gets. Besides the target (modeled as a simple power law) and back¬ 
grounds, all the 2FGL point sources in the field were included in 
the model. 


Most of the sources were undetected on the whole 0.1-100 
GeV energy band. For RBS 1510 (TS=23.0) and 1ES0927+500 
(TS=24.5) we had marginal detections while we could confirm 
RBS 0723, already present in both the 1FGL and 2FGL catalogs. 
Then we computed fluxes in the 1-10 GeV and 10-100 GeV bands, 
relaxing the limit for a detection to TS>9 for the sources that were 
detected in the full band; however, none of the others reached this 
threshold in either energy bin. In absence of a measured flux we 
computed 2 <j upper limits, following iRolke et alj d2005t) . Results 
are collected in Tab. [4] and plotted in Figs. HE In the last two 
columns we also report the number of photons with energy E> 10 
GeV detected within a 0.4° radius (roughly corresponding to 68% 
containment for E > 10 GeV) from the nominal position of the 
source, and the energy of the most energetic one EA rigorous 
study of the significance of these photons, taking into account the 
different PSF of LAT for front and back converted photons, aim¬ 
ing to check and exclude contamination from nearby hard sources 
(though unlikely) and eventually to assess the probability of enclos¬ 
ing background photons within the same aperture, was beyond the 
scope of this work. 

3.4 GALEX data 

The Galaxy Evolution Explorer I Mart in et al]|2005 ) was a NASA 
Small Explorer, in flight since 28 April 2003, and operational until 
mid 2013. 

It performed an all-sky survey in the far UV (FUV, A e fr ~ 154 
nm) and near UV (NUV, A e ft ~ 232 nm) band. We retrieved from 
the MASlJfl archival fluxes observed for the sources of our sample 
(see Tab. [5]l. No data is found for RBS 1176, likely due to a gap in 
the sky coverage of the survey. 


4 SPECTRAL ENERGY DISTRIBUTIONS 

The SED of the sources are reported in Figs. HE In all cases we 
use the same color code: green symbols are used for archival data 
from ASDC, red points for Swift/XKT, UVOT and LAT data and 
black triangles for GALEX data (taken from the database). Note that 

5 Multimission Archive at the Space Telescope Science Institute. 

http://galex.stsci.edu 
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Source Name 

Obs. ID 

V 

b 

u 

wl 

m2 

w2 

BZB J0022+0006 

38113001 

19.23 ±0.17 

- 

19.17 ±0.10 

- 

_ 

19.38 ±0.09 

RBS 0723 

37396001 

17.60 ±0.11 

18.49 ±0.10 

17.65 ±0.08 

17.39 ±0.07 

17.49 ±0.08 

17.43 ±0.05 

1ES 0927+500 

39154001-3 

- 

- 

17.32 ±0.03 

- 

- 

17.97 ±0.04 

RBS 0921 

37547001 

18.78 ±0.33 

19.02 ±0.17 

18.63 ±0.16 

18.79 ±0.12 

18.95 ±0.12 

18.92 ±0.07 

RBS 0923 

48001001-4 

18.97 ±0.19 

19.67 ±0.15 

18.82 ±0.09 

18.59 ±0.09 

18.68 ±0.1 

18.6 ±0.08 

RBS 1029 

36813001-2 

- 

- 

18.47 ±0.06 

- 

- 

18.73 ±0.05 

RBS 1176 

48000001-3 

19.49 ±0.21 

20.05 ±0.15 

19.53 ±0.13 

19.35 ±0.13 

19.10 ±0.13 

19.52 ±0.14 

RBS 1510 

91101001-5 

- 

- 

- 

17.54 ± 0.04 

17.51 ±0.05 

17.65 ±0.04 

RBS 1555 

38303001 

19.37 ±0.11 

- 

19.69 ±0.13 

" 

" 

19.99 ±0.11 


Table 3. Swift/UVOT observed magnitudes for all the archival observations of our sample available until the end of August 2012. Statistical uncertainties only 
are reported: systematic error is always within 0.03 mag and almost generally dominated by statistical ones. 


Source Name 

Fo.l-100GeV 
[10~ 10 ph cm -2 s -1 ] 

r 

TS 

Fl-10GeV 

[10— 10 ph cm -2 s -1 ] 

Fl0-100GeV 
[10 -10 ph cm -2 s -1 ] 

Hard photons 

[Ne> 10GeV] 

Ematr 

[GeV] 

BZB J0022+0006 

- 

_ 

1.88 

< 1.7 

< 0.4 

1 

31.5 

RBS 0723 

11.3±2.6 

1.46±0.08 

79 

3.9± 1.0 

0.69±0.25 

10 

77.6 

1ES 0927+500 

4.2±3.4 

1.38±0.3 

24.5 

< 1.4 

0.62±0.26 

9 

159.1 

RBS 0921 

- 

- 

4 

< 2.4 

< 0.44 

2 

23.4 

RBS 0923 

- 

- 

3 

< 1.5 

< 0.36 

1 

53.7 

RBS 1029 

- 

- 

0.002 

< 2.0 

< 0.29 

0 

- 

RBS 1176 

- 

- 

-0.04 

< 2.3 

< 0.31 

1 

13.7 

RBS 1510 

5.8±3.6 

1.48±0.21 

23.0 

2.0±0.8 

0.36±0.20 

5 

42.9 

RBS 1555 

- 

- 

0.44 

< 2.0 

< 0.39 

2 

10.9 


Table 4. Results of our analysis of the Fermi/LAT data collected until 21 August 2012 for the sources of our EHBL sample. For each target in the first three 
columns the integral flux above 0.1 GeV, the slope of the simple power law model and the TS are reported. Then the value of the flux or the 2<r U.L. in the 
110 GeV and 10-100 GeV energy bins are reported. In the last two columns, the number of photons with E > 10 GeV observed within 0.4° from the nominal 
position of the source, and the energy E maa - the most energetic one. 


Source Name 

FUV Flux 
[AtJy] 

NUV Flux 
[/By] 

BZB J0022+0006 

17.0±0.2 

23.40±0.25 

RBS 0723 

104.3±4.9 

177.2±4.3 

1ES 0927+500 

59.3±5.6 

68.8±4.1 

RBS 0921 

15.1± 1.5 

25.4±1.1 

RBS 0923 

12.6±3.7 

22.8±3.5 

RBS 1029 

17.5±4.3 

32.9±4.9 

RBS 1176 

- 

- 

RBS 1510 

24.9±6.2 

48.1 ±5.3 

RBS 1555 

9.62±0.51 

11.3±0.2 


Table 5. Archival GALEX fluxes in the FUV and NUV bands for our sam¬ 
ple of EHBL. No data is found for RBS 1176, likely due to a gap in the sky 
coverage of the survey. 


in almost all cases the GALEX and UVOT data in the UV filters 
perfectly agree, in spite of the uncorrelated observing epochs. The 
only exception is RBS 1510 for which GALEX provides fainter 
fluxes than UVOT. In this case variability is likely an explanation 
for this difference. 

For comparison, we report in background (gray) the data cor¬ 
responding to 1ES 0229+200 (for simplicity only the observed TeV 
spectrum is shown). It is clear that the structure of the synchrotron 
part of the SED of all the sources closely resembles that of 1ES 
0229+200. In all cases the steep optical continuum is dominated by 
the host galaxy emission. A feature shared by all the sources is the 


large ratio betw een the flux in the UV band and soft X-ray band. 
As re marked in lTavecchio et al J (2009) (see also iKaufmann et al.1 
120111) , in the framework of the one-zone leptonic model, this fea¬ 
tures can be reproduced if the energy distribution of the emitting 
electrons is truncated below a relatively large value, 7 m in = 10 4 ~ 5 . 
As for 1ES 0229+200, this also causes a very hard SSC component, 
consistent with the extremely low flux in the LAT band. 

As discussed in §3.1, for RBS 1176 the observed deficit of 
soft X-ray photons in the XRT spectrum can be interpreted either 
as due to intrinsic or intervening absorption or as the evidence for 
an intrinsic curvature of the spectrum. In the latter case, the ex¬ 
tremely hard soft X-ray continuum (photon index close to 1) would 
be readily interpreted in the SSC scheme as the synchrotron low- 
energy tail of the electrons with Lorentz factor 7 m j n . In this case, 
RBS 1176 would thus be a source characterized by a peculiarly 
large 7 m in, for which the break in the synchrotron continuum - 
which in the other sources occurs between the UV and the X-ray 
bands - is located around 1 keV. Therefore RBS 1176 could be the 
first example of ultra extreme HBL. 

Our selection criterion is rather efficient in selecting EHBL 
blazars, as shown by the similarity of the synchrotron SEDs of the 
selected sources, even if the paucity of the 7-ray data does not allow 
(yet) a detailed comparison of the high energy SEDs. We remark 
that the selection procedure does not explicitly require that the syn¬ 
chrotron peak is in the X-ray band (which is the defining feature of 
EHBL). Nevertheless, almost all the selected BL Lacs show hard 
X-ray spectra, often showing photon indices below 2, locating the 
synchrotron peak energy above 10 keV; it is worth noticing that 
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with indices around 2 the position of the synchrotron peak can shift 
a lot (from ~ 10 16 to ~ 10 19 Hz ), even if the synchrotron lumi¬ 
nosity is almost constant. As shown below, the exact location of 
the synchrotron peak does not affect much the predicted TeV flux, 
since this depends mostly on the low energy edge of the electron 
distribution (see the case of 1ES 0927+500 below). 

The high energy component is unconstrained in almost all 
the sources, with the exception of RBS 0723, RBS 1510 and 1ES 
0927+500, detected by LAT at the highest energies. Based on the 
striking similarity of the synchrotron bump, it is reasonable to as¬ 
sume that the high energy component is also similar to that of 
1ES 0229+200, characterized by an extremely hard s pectru m in the 
LAT (T = 1.36 ± 0.25 according to IVovk et alJ 1201 j. but see 
ICerruti et alj201fj where a softer slope is found) and TeV bands. 
This assumption is supported by the detection of 1ES 0927+500 
and RBS 1510 only at high energy and by the hard photon index 
of RBS 0723 reported in the 2FGL, T = 1.48 ± 0.16. In the next 
paragraph we then estimate the expected TeV flux assuming as SED 
template and physical parameters those of 1ES 0229+200. 


4.1 SED Modeling and predicted TeV fluxes 

We use the one- zone leptonic model of Tavecchio et al.l J 199 8l) . 
fully described in lMaraschi & Tavecchio! 2003). The emitting re¬ 
gion is assumed to be spherical with radius R, filled by a tangled 
magnetic field of intensity B. To reduce the number of free pa¬ 
rameters we assume that the relativistic electrons follow a sim¬ 
ple power law distribution IV (7) = K 7“™ with 7mm < 7 < 
7max- As for th e case of 1ES 0229+200 jTavecchio et alJl2009l : 
iKaufmann et alJl201 ll) this distribution suitably describes the ob¬ 
served SED. The relativistic beaming of the synchrotron and SSC 
radiation is described by the relativistic Doppler factor 8. These 
parameters fully specify the model, that can be uniquely fixed once 
the quantities specifying the two bumps (peak frequencies and lu- 
minosities, spectral slo pes) and the variability timescale are known 
jTavecchio et alJ 19981) . 

In the present case, since in the majority of cases we do not 
have any direct measurement of the high energy hump, it is not 
possible to uniquely derive the physical parameters of the emit¬ 
ting sources. Moreover the minimum variability timescale t var for 
these sources is not known, thus allowing additional degeneracy in 
the K — R plane. On the other hand, relying on the physical pa¬ 
rameters inferred for known EHBL (e.g. jTavecchio et al. 11201 II) . it 
is possible to derive the SSC component and thus an estimate of 
the expected TeV flux. Acco rding to the results of the models in 
ITavecchio et al.l j2010Ll20lIh we then fix the radius R = 6 x 10 15 
cm and the Doppler factor 8 = 20. Then we consider two values 
of the magnetic field bracketing the expected range of the magnetic 
field intensity, B = 0.01 — 0.1 G and for each value we derive 
the remaining parameters (n, K, 7 m in and 7 max ) reproducing the 
synchrotron component (described by the X-ray spectrum and the 
UV data) and the upper limits in the GeV band. The only parameter 
loosely constrained in this procedure is 7max, which, however, has 
only a minor impact on the derived SSC component due to Klein- 
Nishina suppression. For 1ES 0927+500, RBS 0723 and RBS 1510, 
the LAT data allow us to constrain also the level of the SSC com¬ 
ponent and to determine all the parameters. Therefore in this case 
we present only one model. For RBS 1176, as discussed above, we 
consider both possibilities for the X-ray spectrum. For the case as¬ 
suming additional absorption we report two models, as for the other 
sources. For the case in which the X-ray continuum as an intrinsic 
curvature we show one model (black lines), assuming B = 0.01 


G. We do not consider the case of a larger magnetic field for which 

the SSC luminosity (already low) would be m uch smaller. _ 

_ As discussed in lKatarzvnski e t al. (2005) and Tavecchio et all 

(2009), the peculiar SED of EHBL can be reproduced in the frame¬ 
work of the standard one-zone leptonic model assuming that the 
emitting relativistic electrons follow an energy distribution trun¬ 
cated below a relatively large energy or, equivalently, Lorentz factor 
7min- In this case, below the typical synchrotron frequency of elec¬ 
trons with Lorentz factor 7 m in, ^min — 2.8 x Itf Bf^ in 8, the re¬ 
sulting spectrum is described by the characteristics hard power law 
F(u) oc i/ 1 / 3 . The same hard spectrum describes the SSC emis¬ 
sion up to the peak energy hue — 7min m e c 2 8. Since the UV and 
X-ray data constrain Vmin around 10 16 Hz, the typical minimum 
Lorentz factor is 7 m i n ~ 2 x 10 s uff ie BZf 2 S^ 1 ^ 2 , implying 

a SSC peak at uc — 1 u min 16 ^- 2 ^ 8f 2 TeV. Note that in this 
scheme the sources with very large separation between the X-ray 
and the UV fluxes, implying a large value of u m i n , are expected 
to have the maximum of the SSC component at very high energy 
and are thus the most promising sources for TeV detection. Larger 
values of 8 have been sometimes required in the past to model the 
SED of HBL, (e.g. lAleksic et al.l2012h . and would boost even more 
the detectability of our sources. 


4.2 Results 

The resulting theoretical SED are shown by the red (B = 0.1 
G) and blue (B = 0.01 G) lines in Figs. 13151 Solid lines report 
the intrinsic emission, dashed lines show the observed emission, 
corrected for absorption through inter action with the EBL using 
the model of [Dominguez et all j20l ll) . We note here that the ra¬ 
dio emission can not be explained by our model, providing fluxes 
well below the measured level. This is a general feature of single¬ 
region models, which are especially intended to model the emission 
at higher frequencies emitted by compact components. The low fre¬ 
quency radio emission is instead likely produced in extended re¬ 
gions in the jet. 

The cases with low B are characterized by a larger SSC flux 
and viceversa. This is simply due to the well known fact that the 
ratio between the SSC and synchrotron luminosities is proportional 
to the radiation and magnetic energy densities, i.e. Lssc / L syn = 
Urad/UB. For a fixed L syn also (7 ra d is constant, thus Lssc oc 
B~ 2 . 

Figs. l3l5l also display the sensitivity c urves for MAGIC (light 
blue. ISitarek et alj|2013h and CTA ( violet. lActis et al1l201 lh cor¬ 
responding to 50 hours of observation and 5<r significance. It’s 
noteworthy that these sensitivity curves assume a 0.2 dex energy 
binning, therefore dimmer flux densities are still within reach if a 
more coarse binning is adopted. In the low magnetic field case, the 
majority of the sources could be already detected by the present 
generation of_IACTs and all the sources could be easily detected 
by CTA jAcharva et al.ll2013h . In the high B-field case, instead, 
the selected EHBL could be hard to detect even by CTA. How¬ 
ever, we remark that the prediction of the SSC flux for the case of 
high magnetic field, B = 0.1 G, should be considered rather pes¬ 
simistic, since the magnetic field intensity derived for most of the 
known EHBL tends to lie close to the low value B = 10~ 2 G (e.g. 
ITavecchio et al.|[2oTTh . 

For RBS 1176, the case of an intrinsic X-ray break results in 
a large peak frequency (which we recall is directly related to 7 m in, 
that in this case is larger than for the other sources) and a quite low 
SSC luminosity (related to the small energy density of the target 
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Figure 3. Spectral Energy Distribution of the selected extreme BL Lac objects. Green symbols reports historical data from ASI/ASDC archive. Red symbols 
show Swift/UVOT, XRT and Fermi/LXT data discussed in the text. Black symbols display Galex data. Background grey symbols show the SED of 1ES 
0229+200 for comparison. For sources not detected by LAT we report two different models of the SED, corresponding to low (B = 0.01 G, blue lines) and 
high (B = 0.1 G, red lines) magnetic field. For sources with a LAT detection we only report one model, in black (parameters are reported in Tab.[6). For 1ES 
0927+500 w e report two mode ls, corresponding to the two X-ray slopes (see Tab. 2). Dashed lines show the model after absorption with the EBL, calculated 
according to lDominguez et all 12011 ). Light blue and violet curves report the differential sensitivities (5 a, 50 hours of exposure, 5 bins per energy decade) of 
MAGIC and CTA respectively. 


optical-IR photons). Therefore, if the observed lack of soft photons 
is really connected to the intrinsic spectrum RBS 1176 (making 
this source the first ultra-extreme HBL) we do not expect that it can 
easily detected by the CTA. 


5 DISCUSSION 

Extreme HBLs are particular interesting for several reasons, rang¬ 
ing from the study of jet physics to the estimate of EBL and IGMF. 
In this work we have assembled a well defined group of EHBL and 
we have characterized their SED using historical, Swift, Fermi /LAT 
and (when available) GALEX data. All the SEDs closely resemble 
in shape that of the “prototypical” EHBL 1ES 0229+200, demon- 


















Identifying extreme BL Lacs 


9 



46 


45 y 

w 

CiO 

44 

a 

43 cio 
o 


42 


41 



45 

44 

43 

42 

41 

40 


taO 

<D 


& 

tao 

O 



Log v [Hz] 


46 

45 

44 

43 

42 

41 


m 

ao 

u 

a> 


a 

-j 

M 

O 

-J 



a 

J 

DC 

o 


Figure 4. As in Fig. [3] In the case of RBS 1176. the black line refers to the model assuming the highest value of 7 compatible with the X-ray spectrum 
(7 min = 4 X 10 5 , see Tab.[ 6 ]. 


strating the efficiency of our selection in identify these kind of 
sources. We applied a simple one-zone synchrotron SSC model, 
trying to predict the flux at very high energies, showing that most 
of the sources could be detectable by the upcoming CTA. Given 
the shape of the SED, the most effective instruments of the CTA 
to detect EHBL would be the array of Medium Sized Telescopes 
(MST), dominating the overall sensitivity from several hundreds of 
GeV up to few TeV. In case of detection a deep follow up with the 
Large Sized Telescopes (LST) sensitive down to few tens of GeV 
could measure the reprocessed flux and hence the IGMF under the 
hypothesis that the intrinsic emission does not interfere below 100 


GeV. Possible different variability at high (dominated by the intrin¬ 
sic) and low energies could help to separate the two components. 

One of the most interesting features of EHBL is their excep¬ 
tionally hard GeV-TeV spectrum, which in most of cases puts the 
ma ximum of the SSC comp one nt above few TeV. Her e, follow¬ 
ing [KatarzyhskTetal] ( 120051) and lTavecchio et all (1 20091) we inter¬ 
pret the hard TeV spectrum, together with the extreme X-ray/UV 
flux ratio, as the result of the emission from an electron population 
with a large min imum Lorentz factor, 7 m i n ~ 10 5 . As discussed in 
iTavecchio et alj j200gh the exceptionally large 7 m in seem to chal¬ 
lenge existing models for particle acceleration, which usually pre- 
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Source Name 

B 

[G] 

K 

7min 

7max 

n 

BZB J0022+0006 

0.1 

1.7 x 10 11 

3 X 10 4 

2 x 10 6 

3.5 


0.01 

2.75 x 10 13 

9 X 10 4 

6 x 10 6 

3.5 

RBS 0723 

0.15 

6 x 10 8 

2.1 x 10 4 

1.5 x 10 6 

3.0 

1ES 0927+500 1 1 

0.05 

1.7 x 10 10 

4.1 x 10 4 

3 x 10 6 

3.3 

2+3 

0.035 

1.3 x 10 7 

2.7 x 10 4 

3 x 10 6 

2.7 

RBS 0921 

0.1 

8 x 10 7 

4.7 x 10 4 

1.8 x 10 6 

2.8 


0.01 

6 x 10 9 

1.3 x 10 s 

6 x 10 6 

2.8 

RBS 0923 

0.1 

1.4 x 10 8 

2.2 x 10 4 

2 x 10 6 

3.3 


0.01 

5.1 x 10 12 

1.2 x 10 s 

5.2 x 10 6 

3.3 

RBS 1029 

0.1 

1.4 x 10 8 

2.2 x 10 4 

2 x 10® 

3.0 


0.01 

1.4 x 10 10 

7 x 10 4 

6 x 10® 

3.0 

RBS 1176 2 

0.1 

6 x 10 7 

2.3 x 10 4 

10® 

2.8 


0.01 

4.6 x 10 9 

7.3 x 10 4 

3 x 10® 

2.8 


0.01 

3 x 10 11 

4 x 10 5 

3 x 10® 

3.1 

RBS 1510 3 

0.12 

6.2 x 10 9 

2 x 10 4 

2 x 10® 

3.35 

RBS 1555 

0.1 

1.2 x 10 6 

1.3 x 10 4 

3 x 10® 

2.6 


0.01 

7.5 x 10 7 

4.3 x 10 4 

10 7 

2.6 


Table 6 . Physical parameters describing the jet emission in our sample whithin the framework of the SSC model from iTavecchio etal (1998), with the 
simplifying assumption that the energy spectrum of the emitting electrons follows a simple (instead of broken) power law distribution. A common Doppler 
factor 5 = 20 is assumed for all the sources. For each target source the intensity of the magnetic field within the emission region B, the electron spatial density 
K, the minimum and maximum electron energy and 7 max, the slope of the power law distribution n are reported. 1 : R = 6.5 X 10 15 cm, S = 30; 

2 the first two lines refer to the model assuming a power law X-ray spectrum. The third line is for a broken-power law X-ray spectrum (black line in Fig. 4); 
3 :R = 1.4 x 10 16 cm. 
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Figure 5. As Fig.[3] 


diets power-law non-thermal ta ils starting from much lower values 
(e.g. lSironi & Spitkovskyfeoi lh . 

Another intriguing point is the apparent stabili ty of the TeV 
spectrum over s everal years jAharonian et al][20071 . but see also 
lAliu et al1l2014l for the case of 1ES 0229+200 where marginal ev¬ 
idence for mild TeV variability on a several months/annual time 
scale is found), which neatly distinguishes EHBL from HBL that 
show rather variable TeV emission, down to timescales of few min¬ 


utes. As discussed in ITavecchio et al] (120091 ) for 1ES 0229+200 
the radiative cooling time of electrons at 7 m in is of the order of 2 
years, barely consistent with the observations. However, it is dif¬ 
ficult to avoid other sources of losses, such as adiabatic losses. 
One interesting possibility to understand such a stable spectr um 
(e.g. iMurase et al.ll20l 2al : lEssev et alj|201 ll : iTavlor et aPl201 ill is 
that the GeV-TeV continuum is not the primary emission from the 
source, but instead is reprocessed radiation from cascading pro¬ 
cesses spread into the intergalactic space. Electromagnetic cas¬ 
cades could be the result of the absorption of primary multi-TeV 
gamma rays or, instead, could be produced by Bethe-Heitler pair 
creation or photo-meson reactions involving ultra-high energy pro¬ 
tons accele rated in the source and beamed toward the Earth. As 
shown in ITavecchio! d2014l) . the physical parameters of the jets 
are consistent with the requests of the hadronic cascade scenario, 
both in te rms of maximum had ron energy and je t powe r. As dis¬ 
cussed by IMurase et al.l d2012al) and lTakami et ~al] d2013l) an effec¬ 
tive test to distinguish between intrinsic and reprocessed emission 
is the observation of photons at several TeV, only possible in the 
latter case, since the cascade emission, being produced at lower 
distance, is less affected by absorption than the intrinsic one; al¬ 
though e xotic process es such as the photon-axion conversion (e.g. 
Ide Angelis et alj201 lb could result in a lower effective absorption. 
Due to absorption and reprocessing of primary TeV photons a sig¬ 
nificant contribution to the GeV extragalactig gamma-ra y back ¬ 
ground (EGB) could be re l ated to E-HBL (se e e.g. Ilnoue & Iokal 
120121 : IMurase et al1l2012bl ; lAiello et aT]|2014l) . We plan to study 
this issue in a forthcoming paper. 

Another topic worth of discussion is the position of these 
EHBLs within the general BL Lac population and the nature of 
their parent population. By construction, the sources we are dis¬ 
cussing have very faint radio emission (radio luminosities around 
L r ~ 10 4u erg s~ 4 ) indicating very weak radio jets (e.g. 

IGiroletti et al] 120041 . for some sources). Considering that at least 
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a fraction of the radio flux comes from the beamed jet synchrotron 
component, we argue that the intrinsic radio luminosity drops be¬ 
low the lower end of the FRI radiogalaxy power range, making 
EHBLs suitable candidates for the aligned counterpart s of th e weak 
radiogalaxies population studied bv IBaldi & Capettil (2003). 

Finally we wish to speculate on the existence of even more 
extreme EHBL. In fact it is tempting to suppose the existence of 
blazars which SED is shifted toward even higher energies, with the 
X-ray emission peak above 10 keV and high-energy component 
reaching several tens of TeV. Given the very high energy of the 
synchrotron and IC peaks and low flux levels these sources could 
escape detection by Swift and Fermi. We speculate that possibly 
the CTA survey will allow to detect or constrain this hypotetical 
subclass of the HBL population in the TeV band. Of course such 
sources would deeply suffer from EBL absorption around their IC 
emission peak, limiting their detectability at low redshifts and the 
detection of the left side of the IC bump below few TeV. Arguably, 
a relatively strong fraction of their radiated power would be repro¬ 
cessed into the GeV band, setting limits to their spatial density. We 
will further investigate this issue in a future work. 

A generalization of our criterion, for instance relaxing the con¬ 
dition of the known redshift, will allow to validate this EHBL se¬ 
lection method on new and larger samples. This class of sources 
is especially important for improving our understanding about the 
far-IR EBL, IGMF and EGB. Exploitin g at th e same time the ca¬ 
pabili ty of NuStar l Harrison et al] 2013b SKA JCarilli & Rawling si 
l2004h and CTA dAcharva et alJ 2013b will be particularly reveal¬ 
ing. Even b efore the comple tion o f CTA, the plann ed ASTRI/CTA 
mini-array jDi Pierro et al]|2ol3l : 1 Vercellone et alj|2013h could be 
exploited in this direction. Actually, whereas the brightest flux from 
this sources is expected to be in the band up to 1 TeV, where MSTs 
will be the most sensitive instruments of CTA, still crucial physical 
information is engraved in the high energy tail of the IC bump, in 
the band best observed with the Small Sized Telescopes. 

For instance, these extreme sources could allow eventually to 
compare leptonic and hadronic emission scenarios in a multi-TeV 
territory where degeneracy of competing models is significantly re¬ 
duced (see e.g. iMurase et al.ll2012ab . Opportunity would arise for 
tests of non standard physics such as violations of the Lorentz 
invariance (Fairbairn et al] |2 0 14h or the_speculated existence of 
axion-like particles ( Meyer & Conradll2014b . once the number of 
these multi-TeV photon factories should become greater than now. 
Therefore EHBL will be very interesting targets for challenging hot 
topics in fundamental physics, both with the ASTRI/CTA Mini- 
Array (improving the current H.E.S.S. sensitivity above 5 TeV) and 
with the full array of CTA small telescopes, which dominates the 
overall CTA sensitivity above the same threshold. 
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